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We consider the problem of option pricing under stochastic volatility models, focusing on 
the linear approximation of the two processes known as exponential Ornstcin-Uhlcnbcck 
and Stein-Stein. Indeed, we show they admit the same limit dynamics in the regime of 
low fluctuations of the volatility process, under which we derive the exact expression 
of the characteristic function associated to the risk neutral probability density. This 
expression allows us to compute option prices exploiting a formula derived by Lewis and 
Lipton. We analyze in detail the case of Plain Vanilla calls, being liquid instruments for 
which reliable implied volatility surfaces are available. We also compute the analytical 
expressions of the first four cumulants, that are crucial to implement a simple two steps 
calibration procedure. It has been tested against a data set of options traded on the 
Milan Stock Exchange. The data analysis that we present reveals a good fit with the 
market implied surfaces and corroborates the accuracy of the linear approximation. 

Keywords: Econophysics; Stochastic Volatility; Monte Carlo Simulation; Option Pricing; 
Model Calibration 
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1. Introduction 

The recent financial crisis has emphasized the need for reliable quantitative analysis 
of market data, able to guide the formulation of realistic theoretical models for the 
dynamics of the traded assets. The Black-Scholes (B&S) and Merton approach 
to option pricing |l|2j assumes a Gaussian dynamics for the underlying assets and 
therefore it fails to reproduce the well known stylized facts exhibiting clear evidence 
of deviations from the normality assumption. This is why, in recent years, more 
realistic alternative models have been proposed in the literature. In particular, to 
capture the time varying nature of the volatility, assumed to be constant in the B&S 
approach, the stochastic volatility models (SVMs) represent a theoretical framework 
used both in the research and the financial practice. Among the most popular SVMs 
include the Heston gj, Stein-Stein (S2) @], Schoblc-Zhu [5], Hull- White [B] and 
Scott [7] models. For reviews of SVMs we refer to |8I9I10| . More recently, the model 
known in the econophysics literature as exponential Ornstcin-Uhlcnbcck (ExpOU) 
has drawn particular attention because of its ability to reproduce a log-normal 
distribution for the volatility, the so called leverage effect as well as the evidence of 
multiple time scales in the decaying of the volatility auto-correlation function [llj . 
In [12] , the statistical characterization of the process under the objective probability 
measure has been carried out from both the analytical and numerical points of view. 

As far as the pricing problem is concerned, semi closed-form expressions for 
the price of European options are available for the Heston and the S2 models. For 
the ExpOU model the problem was addressed in the original paper by Scott, who 
worked out a quasi closed- form pricing formula in the spirit of [6]. However, it even- 
tually relies on the Monte Carlo (MC) simulation of the history of the volatility, 
while prices under Heston and S2 can be efficiently computed exploiting Fast Fourier 
Transform (FFT) numerical techniques |13|14|15|16j . Recently, based on the Edge- 
worth expansion of the risk neutral density, an analytical expression of the pricing 
function for ExpOU has been derived |17|18j . The accuracy of their approach has 
been tested numerically, and, at least for the considered regime, the approximate 
probability density function (PDF) they provide is unable to fit the one recon- 
structed via MC [12]. In conclusion, a satisfactory solution to the pricing problem 
under the ExpOU model is still lacking and it deserves further investigation. Indeed, 
the aim of the present paper is to discuss, under a risk neutral framework, its linear 
approximation for which a complete analytical characterization of the character- 
istic function (CF) can be provided. This allows to employ the Lewis and Lipton 
formula to efficiently solve for derivatives prices and to test the accuracy of the 
approximation in reproducing market observed volatility smiles. 

The paper is organized as follows. In Section[2]we review the risk neutral formu- 
lation of a SVM when the dynamics of the stochastic variable driving the volatility 
is described by the Ornstein-Uhlenbeck process, as for the ExpOU and the S2 mod- 
els. We show how, under the regime of low fluctuations of the volatility process, the 
returns dynamics reduces to a linear one and we derive the exact analytical expres- 
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sion of the corresponding CF. As presented in | Appendix A the explicit expressions 
of the first four cumulants, whose knowledge allows for an efficient calibration pro- 
cedure, has been computed. In Section [3] we perform a cross-sectional fitting of the 
Linear, ExpOU and S2 models, evaluating the parameters from a data set of Plain 
Vanilla call options and detailing the steps of the adopted calibration methodol- 
ogy. The ability of the Linear model to reproduce the data and the accuracy of 
the approximation are evaluated comparing the volatility smiles reconstructed after 
calibration with the original market ones and with those exhibited by the ExpOU 
and the S2 models. The final Section draws the relevant conclusions and suggests 
some possible applications of the analytical results in the field of market risk man- 
agement, such as Valuc-at-Risk and Expected Shortfall evaluation. 



2. The Linear Model 

The class of SVMs we consider is described by the following system of stochastic 
differential equations (SDEs) 

dS(t) = pS(t)dt + a(Y,t)S(t)dWt(t) , S(t ) = S Q ; 

dY(t) = a(7 - Y{t))dt + kpdW x {t) + k^l-p 2 dW 2 {t) , Y(t Q ) = Y , (2.1) 

where dWi and dW 2 are two independent Wiener processes, while So, fi, Yq, a, 7, k 
and p are constant parameters. The dynamics of Y(t) corresponds to an Ornstein- 
Uhlenbeck process, whose stationary mean and variance are given by 7 and (3 = 
k 2 /(2a). The rate of convergence to the steady state is given by 1/a, while the 
correlation parameter p takes value in [—1,1]. The volatility a(Y, t) is a smooth 
function of Y and t and defining a(Y, t) = me y ^ we obtain the ExpOU model 
(TUT], while for a(Y,t) = mY(t) the S2 model gE] is recovered. 

Given the market model specified by Eq. ()2.1|) . the standard approach to option 
pricing consists of passing to an equivalent risk neutral measure P* under which 
the discounted price process S{t) = e~ rt S(t), with r the risk-free interest rate, is 
a martingale. Indicating with E*[-] the expected value under P*, this martingale 
property simply reads E* [5(t)|S'(to)] = S(to) and the risk neutral dynamics of the 
model becomes 

dS(t) = rS(t)dt + a(Y,t)S(t)dW?(t), 

dY(t) = [a(7 - Y(t)) - kr)(S, Y,t)] dt + kpdW*(t) + ky/l - p 2 dW;(t). (2.2) 

In the above equation, W* (t) and IU 2 * (t) are independent standard Brownian mo- 
tions under the measure P* and the function r/(S,Y,t) correcting the drift term of 
Y(t) is called the market price of volatility risk, see [19] . The function r\ depends on 
the variables t, S, Y and not on the contract parameters. It takes the same form for 
different derivative contracts stipulated on the same underlying S(t), parametrizes 
the space of risk neutral measures and defines a(j — Y) — krj{S,Y,t) which is the 
risk neutral drift of Y . 
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S and Y being Markovian processes, rj is a function of the processes at time t,rj = 
r](t,S(t),Y(t)); apart from suitable intcgrability conditions, from a mathematical 
point of view rj is an arbitrary function and we assume it to be a linear function of 
the process Y(t) 

V (Y)= Vo +r h Y(t). 

In the light of its arbitrariness, our choice of a linear rj is eventually dictated by the 
opportunity to preserve the mean reverting Ornstcin-Uhlcnbeck dynamics for Y in 
the following, and not by any financial intuition. It is worth noticing that this choice 
applies to the whole class of models (|2.1|) . independently on the explicit functional 
form of a(Y,t). It is coherent with the assumption done in |18j and generalizes the 
one made by Stein and Stein in their original work [4] corresponding to r/i = 0. We 
can redefine the parameters a and 7 as 

, ; a aj-kr/o 

a — > a = a + kr]i , and 7 — > 7 = , 

a 

and, after applying Ito's Lemma to the centred logarithmic return X(t) = In S(t) — 
lnS'(to) — r(t — to): finally the risk neutral dynamics reads 

dX(t) = --a 2 (Y,t)dt + a(Y,t)dW*(t); 

dY(t) = 5(7 - Y(t))dt + kpdW*(t) + ky/l -p 2 dW;(t), (2.3) 

with initial conditions X(to) = and Y(to) = Yq, and a > that ensures the 
stationarity of the Y process. 

When the stationary variance of Y is small, /3 = k 2 /(2a) <C 1, we can perform 
a first order Taylor expansion of a and a 2 around Y — 7. Defining the process 
Z = Y + 1 — 7, and the parameters m = me 7 , k = k for the ExpOU model, and 
Z = Y/j, rh = mj, k = fc/7 for the S2, the processes in (|2.3|) reduce to 

~ o 

m 

dX(t) = — —(2Z(t) - l)dt + mZ(t)dWf , X(t ) = X = 0, (2.4) 

dZ(t) = a(l - Z(t))dt + kpdWKt) + kyj\ - p 2 dW£(t) 1 Z(t Q ) = Z . (2.5) 

The accuracy of the approximation has been discussed in detail for the ExpOU 
model in \12\ . In particular, a numerical analysis based on MC simulation for /? < 
10% supports the linearisation leading to previous equations. However, the linear 
approximation docs not preserve the martingality of S(t), which is a crucial point 
for pricing purposes. From the definition S(t) = S(to)e xi - t \ the violation of this 
property is readily assessed by computing the deviation of E*[e x (*)|Xo] from 1. In 
order to obtain a martingale dynamics we modify the drift term driving X(t) by 
means of a deterministic time dependent function A4(t — to). The corrected SDE 
reads 

~ 2 

Tfl 

dX(t) = — — (2Z(t) - 1 + M(t - t ))dt + mZ(t)dW* . (2.6) 



Option pricing under Ornstein- Uhlenbeck stochastic volatility: a linear model 5 



The returns transition probability distribution p x (x, t\Xq, Zq), with r = t — to, can 
be expressed in terms of the CF /(</(>, t;Xo, Zq), implicitly defined by 

e-^/(0,r;X o ,Z o )d0 . (2.7) 



1 f A 

p x {x,r\X ,Z a ) = — J 



The CF satisfies the Fokker-Planck backward equation in the Fourier space associ- 
ated to the two-dimensional process described by Eqs. Q2.50 and (|2.6[) . This equa- 
tion can be solved exactly, following a standard technique (see |3l20j ). by guessing 
a solution of the form 

f m 2 r 

f{4>, r; X , Z ) = cxp j -icp— j X(r')dr' + A{<f>, r) + B(0, r)Z + C(4>, t)Z 2 + i<j>X 

(2.8) 

The explicit expressions for the three functions A, B and C represent the main 
analytical result of the paper and read El 



a(4>,t) 



2a- 



2k 2 



n — h 



b-d 



1 



-dr\ 



2 

2k 2 



[In (l - ge 
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(l+ge-i^ (e-* T -l) 



2 

6-d 1 



(1 - .9) (1 - ge- d -) 
h 



d(l - ge~ dT ) 



-dr 



(2.9) 
(2.10) 
(2.11) 



Ak 2 l-ge- dT ' 

where we have introduced the auxiliary functions b = 2a(l—ip<&), d = \/2a 2 $> 2 + b 2 , 
g = (6 - d)/(6 + d), h = iam^/k. n = a(b - d)/{2k 2 ). and $ = kmd)/a. It is 
relevant noting that the difference between principal logarithms in the second line 
has not been contracted into the logarithm of the ratio. Indeed, this operation can 
be performed only by taking into account a suitable correction (see Eq. (2.4) in [2T|V 
In order to assign the function M., we impose E*[e x O|X ] = f(—i,T;X ,Z ) = 1 
thus finding 

M(t) = 4r-^- [A{~i,t) + B{-i,t)Z q + C{-i,t)Z 2 +X q ] . (2.12) 
m z dr 

Previous expression a posteriori justifies the choice of M. as an homogeneous func- 
tion of time in Eq. (|2.6[) . 



a From now on we shall drop the tilde over the model parameters. 
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3. Numerical Results 

3.1. Cross-sectional fitting 

In the financial practice, SVMs are first calibrated on market data and then used 
for pricing. The calibration of the model parameters could be performed following 
different approaches and this problem has been widely addressed in the literature. 
Several procedures have been proposed in different contexts, see e.g. j9!22l23j and 
a discussion concerning asymptotic formulae for the implied volatility smile can be 
found in |24j . which also provides an overview of asymptotic methods. In this work, 
we exploit the relationship between implied volatility smiles and the variance a^, 
skewness ( T and kurtosis n T of the risk neutral PDF, as provided by the expression 
given in [25] (see also |26l27j ) 



'"imp,r(^l) — / — 



(3.1) 

where <7i mPiT is the B&S implied volatility for the time to maturity r and 



The expression p.l[) is based on the approximation of the risk neutral PDF for 
fixed r by means of a Gram Charlier expansion and it does not rely on the choice 
of any specific underlying dynamics. Its range of applicability is discussed in detail 
in [25] . where it is shown that Eq. (|3 . 1 [) is effective for d\ ~ and a T <C 1, which is 
realized in practice (for r ~ 1 year, typically a T ranges from .2 to .3). 

We can compute the cumulants for the Linear model exploiting the analytical 
formulae provided by Eq. (|2.8[) - (|2.1ip . to which they are related through 



., ra 9"ln/(0,r;Xo,Z o ) 



(3.2) 



4=0 



To obtain the analytical expressions, using MATHEMATICA® we approximate the 
logarithm of / by means of a 4-th order Taylor expansion around = and then we 
extract the four coefficients of the expansion and multiply them by the appropriate 
constant factor, finally finding the results reported in | Appendix A| After identifying 
o\ with fc2,T) then the skewness and kurtosis read £ r = k^^/a^., and n T = k^^ T jo\, 
respectively. The exact analytical CF is also available for the S2 model (see |4|5j ) 
and following the same approach it is possible to compute explicitly the related 
expressions for the cumulants. 

We limit our analysis to Plain Vanilla call options, whose implied volatilities 
Cimp.r are available from market data providers for different maturities r and strike 
prices K. The underlying spot price So and the term structure of risk- free rates 
r r can be retrieved from the market as well. Table 1 sums up the complete data 
set available, corresponding to options written on the Intesa San Paolo S.p.A. asset 
with spot price S = 5.16 EUR, as of 22nd November 2007 on the Milan Stock 
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Table 1. Implied volatilities market data. 
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Exchange. Annualized implied volatilities values are quoted, with the corresponding 
log-moneyness, time to maturities and risk-free rates retrieved from the EUR yield 
curve. The complete list of parameters to calibrate is given by m, Yq (equivalently 
Z ), fc, a, 7, and p, see Eq. (|2.3[) . Under the hypothesis that the process Y(t) driving 
the volatility has reached the stationary state, we fix Yq = 7 with 7 = for the 
Linear and ExpOU models, and 7 = 1 for S2, see also [T2]. In order to fit the 
remaining four free parameters, we adopt the calibration procedure detailed below. 
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Fit o~ T , £ T and k t from market smiles. Exploiting Eq. (|3.1[) . which provides an 
approximation of the smiles in a suitable region around d = 0, we fit the empirical 
data with a Marquard-Lcvenbcrg algorithm [28] retrieving the optimal values a^ 1 k , 
(Jf , and n\ lk . These values and the associated standard errors e^ fe , , and e^ fe , 
are summarized in Table 2. The last column shows that smiles made of fewer points, 
like the one corresponding to r = 0.8274 (4 points), result in a greater estimation 
error, but we expect this effect to be greatly reduced for larger data sets, if available. 
In Fig. Q] we present the market data and the parabolic approximation Eq. (|3.ip 



Table 2. Market calibrated normalized cumulants and their 
standard errors. 
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with parameters fixed as in Table 2. Individual smiles are very well reproduced 




Fig. 1. Implied volatilities for Intesa San Paolo: market data (dark points) and parabolic approx- 
imation, sec Eq. QO) . 



and for long time to maturities the curves flatten, as expected, while the highest 
implied kurtosis corresponds to the shortest r. For r = 0.3260 yr we notice that 
the volatilities for extreme positive log-moneyness are suspiciously out of scale. We 
calibrate on the entire data set, but we expect that this large fluctuations will not 
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be reproduced by the models under investigation. 

Find optimal a, k, m, and p from the time scaling of a T , C, T and k t . The cali- 
bration can be done by fitting the values reported in Table 2 with those computed 
from the models. The scaling of er r , £ r and k t with r is known analytically for the 
Linear and S2 models, while it can be estimated numerically for ExpOU. For the 
latter case, we sample Nmc = 10 5 paths and we obtain the MC estimators cr^ 1 , 
Ct IC > k t IC with associated standard errors e^ c , e^ c and e^ c . The set of optimal 
values satisfies the equation 

(„Mk „ \2 (/-Mk /* \2 I u .Mk \2 

,2 + 2 + 2 

e <7 T £ C T % 

where e 2 ^ = e^ feZ + e^ C ^ an d analogously for e 2 anc j e 2^^ ^ ne jyjQ error being 
zero for the Linear and S2 models since their cumulants are known analytically 
and no MC simulation is required. The optimization problem is solved by means of 
MINUIT routines [53] and the final results are contained in Table 3. The value of (3 
of order 10% from the calibration of the ExpOU dynamics supports the accuracy of 
the linear approximation, while no statistically significant differences are observed 
between the Linear and the S2 models. 



, to* , p* = argmin 

a,k,m>0, (»£[— 1,1] T 



Table 3. Optimal parameters for the three models. 
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0.1 
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0.018 


-0.51 


0.09 


0.13 


0.04 


S2 


5.7 


1.3 


1.9 


0.4 


0.265 


0.008 


-0.41 


0.07 


0.32 


0.14 


Lin 


5.6 


1.3 


1.9 


0.4 


0.264 


0.008 


-0.41 


0.07 


0.34 


0.15 



As a final comment, we point out that the most computationally instcnsive step 
of the above procedure corresponds to the MC simulation of the ExpOU dynamics, 
which was performed on a i686 machine equipped with Intel(R) Core(TM)2 Quad 
CPU Q6600 @ 2.40GHz processor, taking about 1500 minutes to complete. Taking 
advantage of the analytical cumulants, the calibration of the other two models does 
not suffer this limitation and requires only a few seconds. 



3.2. Option prices and implied volatility smiles 

With the parameters values reported in Table 3 we compute the option prices, 
extract the implied volatilities and plot the reconstructed volatility smiles against 
the market ones in order to asses the accuracy of the Linear model. We also compare 
the results with the smiles generated by the ExpOU and S2 models. 

The option prices for the ExpOU model have to be MC computed, while in the 
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other cases the knowledge of the CF allows to use the pricing formula [14115] 

c(s ,t ) = -fv c r +oo dz it^- 



2tt 

St 
2ir 



ic—oo 
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+00 

duj sin(wL>) lm[W(oj+)f(-ui+) - W(u}-)f(-u)-)] 







(3.: 



where we have dropped the dependence on r, Xq and Zq in the characteristic 
function /. In the previous formula, D = log (Sq/ K)+rr, Sq is the spot price, K the 
strike, W(ui) = [uj 2 — iuj]~ l , lu+ = +ui + ic, oj- = —uj + ic, and z = uj + ic. When p S 
( — 1, 1), for the Linear model Eq. (|3.3|) applies for c £ {c : c > l}n{c : c_ < c < c + }, 
where = a/[km(p^f 1)] correspond to the imaginary part of the singularities of 
the CF (|2.8j) . For p ~ —0.5, a ~ 6, k ~ 2, and m ~ 0.26 the above intersection 
is not empty and we set c = A^j-^ with A = 0.5 (we have verified that our 
results are insensitive to different choices of A 6 (— (1 + p), 1)). The identification 
of the singularities and a similar analysis has been performed for the characteristic 
function of the S2 model. 

We compute the integrals involved in Eq. (|3.3[) using an adaptive trapezoidal 
algorithm, optimized to calculate sine or cosine transforms, e.g. see the routine 
dqawf .f available at http://www.netlib.org/quadpack/. In order to check the re- 
liability of the above numerical setup we have compared the implied volatilities 
obtained through inversion of the prices computed via (|3.3p with those from a MC 
simulation (Nmc = 10 7 , parameters fixed as in Table 3) of the Linear model. For 
all the available time to maturities, the implied volatility smiles and the simulated 
ones are in full agreement at 68% confidence level. 

In Fig. [2] we present the implied volatilities smiles reconstructed by the Linear 
model, showing its ability to capture the correct shape of the volatility. The agree- 
ment slightly decreases for deep in and out of the money options; from the third 
panel corresponding to r = 0.326 yr we notice that the fluctuations of the two 
outlying points are not reproduced, as expected. The statistical uncertainty on the 
parameters values reflects in the standard errors associated to the curves. The error 
propagation has been performed making use of the first order derivatives computed 
by means of finite difference methods. It is worth pointing out that, following the 
guidelines depicted in the previous Section, calibration and pricing are carried out 
in an efficient way; this is mainly due to the analytical characterization we provided 
for the CF and cumulants, making the entire approach a real time procedure. 

In Fig. [3] we present the smiles for the set of parameters corresponding to the 
ExpOU calibration and we plot both the curves obtained from MC simulation 
(Nmc = 10 5 ) of the exponential model and those computed integrating Eq. (|3.3[) 
for the Linear model with the same parameters values (curves corresponding to the 
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Fig. 2. Comparison between market implied volatility smiles for Intesa San Paolo and those ob- 
tained from the Linear model. Parameters values as in the third row of Table 3, as obtained from 
the calibration of the Linear model. 



linear case have been shifted rightward). In Fig. [H we plot the PDF of the Linear 
model against the one for the ExpOU model computed with trapezoidal integra- 
tion and MC simulation, respectively. Both panels confirm the analysis performed 
in [T2], showing fatter tails and a lower central peak for the histogram of the ExpOU 
with respect to the PDF of the Linear model. Even though the value f3 ~ 13% (see 
Table 3) is at the edge of the regime allowing the linearisation, as far as the volatility 
smiles obtained from the two models are concerned, we conclude that they are in a 
good statistical agreement. In Fig. [5] a comparison analogous to the one in Fig. [3] 
for the S2 parameters is reported, revealing again the statistical agreement. With 
respect to Fig. [3l the narrower error bars reflect the fact that parameters fitting has 
been performed exploiting the available analytical information. Actually, the MC 
simulation involved both in the calibration and the price computation for ExpOU 
introduces an additional statistical uncertainty. 



4. Conclusions and Perspectives 

This paper deals with the problem of option pricing under stochastic volatility 
and calibration to market smiles. The main focus is on a specific SVM where the 
dynamics of financial log-returns is driven by linear drift and diffusion coefficients, 
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Fig. 3. Comparison between market implied volatility smiles for Intesa San Paolo and those ob- 
tained from the ExpOU and Linear models. Parameters values as in the first row of Table 3, as 
obtained from the calibration of the ExpOU model. 
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Fig. 4. Comparison between p x (x, r\Xo, Zo) with r = 1 for the Linear and ExpOU models in 
log-linear (left panel) and linear (right panel) scales. Parameters values as in the first row of Table 
3. 



which we show to be the limit case of the exponential Ornstein-Uhlenbeck and 
the Stein-Stein models under the low volatility fluctuation regime. The analytical 
contribution we provide is the exact characterization of the characteristic function of 
the Linear model under risk neutrality. Following the market practice, we calibrate 
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Fig. 5. Comparison between market implied volatility smiles for Intesa San Paolo and those ob- 
tained from the S2 and Linear models. Parameters values as in the second row of Table 3, as 
obtained from the calibration of the S2 model. 



the model on a sample of implied volatilities by means of the two step calibration 
procedure detailed in Section [3J In this regard, the knowledge of the analytical 
expressions of the cumulants substantially reduces the computational effort required 
to fix the parameters values. For the considered data, we are able to quantitatively 
asses the capability of the model to reproduce the market volatility smiles, finding a 
statistically significant agreement with the empirical curves. By means of the same 
procedure, we also compute the implied volatility values under the ExpOU and S2 
models after calibration to evaluate the accuracy of the linear approximation. In 
both cases, the propagation of the statistical uncertainty results in an error band 
which clearly reveals that the three models are all in full agreement. In particular, 
the measured value of /? for the ExpOU model justifies its linearisation, whose 
degree of analytical tractability makes it more desirable in view of option pricing. 

In this work, we exploited the Greeks to evaluate the statistical uncertainty of 
the reconstructed volatility smiles; as a future perspective, it would be interesting to 
analyze the sensitivity with respect to movements of market volatility curves, which 
is what is done in financial practice by traders to hedge option positions. We also aim 
at applying the Linear model considered here in the context of market risk measures, 
like Value-at-Risk and Expected Shortfall, exploiting the knowledge of the analytical 
CF following the guidelines traced by [3D]. For risk management purposes, it would 
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be interesting to analytically characterize the PDF tails. Indeed, the numerical 
convergence of the integral in Eq. (|3.3j) excludes a power law decay, but an explicit 
analytical result such as for the Heston [31] and the ExpOU models [32] discerning 
between exponential, Gaussian or different scalings is still lacking. 
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Appendix A. Linear model cumulants 

In the following we report the analytical expressions of the cumulants of the Linear 
model. 
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